Specific heat of classical disordered elastic systems 
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We study the thermodynamics of disordered elastic systems, applied to vortex lattices in the Bragg 
glass phase. Using the replica variational method we compute the specific heat of pinned vortons in 
the classical limit. We find that the contribution of disorder is positive, linear at low temperature, 
and exhibits a maximum. It is found to be important compared to other contributions, e.g. core 
electrons, mean field and non linear elasticity that we evaluate. The contribution of droplets is 
subdominant at weak disorder in d = 3. 

PACS numbers: 



Understanding the temperature dependence of the spe- 
cific heat in glasses remains puzzling, e.g. the linear 
low T behaviour observed in structural glasses |1J and 
spin glasses In some systems the crossover tempera- 
ture from quantum to classical behaviour may be quite 
low. The phenomenological two levels system model Q 
yields a linear behaviour both in the classical and quan- 
tum regime. Although the classical problem appears sim- 
pler there are only few, and mostly mean field, solvable 
models of glassy systems where one can actually compute 
the specific heat |2j, |^ . Ordered systems with continous 
symmetry admit spin wave type excitations which yield a 
T-independent specific heat C v = C eq from the equipar- 
tition of the energy. Non linearities such as quenched 
disorder will cause a deviation which is interesting to 
characterize and compare with the linear contribution 
from the two level system arguments. 

A class of glassy systems recently much studied are 
disordered elastic systems, ranging from vortex lattices 
[EL H, electron crystals ||, charge and spin den- 
sity waves 10] to disordered liquid crystals In all 
of these systems the competition between disorder and 
elasticity leads to pinning and glassy behaviour. Spe- 
cific heat measurements in density waves gapped systems 
showed intermediate linear and sublinear regimes with 
non-equilibration effects [L2I Il3j . In superconductors in 
a field, the H and T dependence of C v relates both to 
the symmetry of the order parameter and to the thermo- 
dynamics of the vortex lattice. If the contribution of the 
normal electrons in the vortex cores dominate, the stan- 
dard expectation is that C v is linear in T with a linear in 
H dependence for s wave [lij and ff 1 / 2 for ci- wave [l5j . 
A specific heat linear in temperature has indeed been 
measured in various materials. A 7? ^-dependence has 
been observed and arg ued for d-wave superconductivity 
in YbaCuO [IrilrRll3 |. However the nonlinear H depen- 
dence observed in other, a priori non <i-wave materials, 
is a well known puzzle, as discussed in [HHi. On the 
other hand, the contribution of the phonons of the vortex 
lattice VL, the so-called "vortons" , seems to be within ex- 
perimental resolution |2l|. and may lead, within the full 



temperature range below VL melting, to more complex 
behaviors. These were analyzed in absence of disorder, 
assuming a dissipative quantum dynamics with friction 
r] arising from interactions between vortons and vortex 
core electrons |22], |23J. It yields again C v oc rjT with 
different H dependence but only for T < , the vor- 
ton Debye-like temperature which is poorly known (esti- 
mates in YBaCuO ran ge fr om well below IK up to 10K 
in the superclean limit |22j). Above one recovers the 
equipartition value to which the specific heat anomaly 
was compared at melting |2J, |25|, 1261 127| . These analysis 
however neglect disorder and other non linearities which 
for T > Tjf can be treated classically. 

In this Letter we compute the specific heat of an elas- 
tic system in presence of pinning disorder in the classical 
regime. We show that disorder produces a substantial 
rise above equipartition, linear at low temperature and 
exhibiting a maximum at a characteristic depinning tem- 
perature. We show that in d > 2 the contribution from 
the two well droplet arguments is subdominant at weak 
disorder. We find that the disorder contribution is quite 
sizable compared to other contributions, e.g. of the non 
linear elasticity that we also evaluate. These results hold 
for a periodic object, i.e. a Bragg glass, or for interfaces 
with continuous degrees of freedom. In a companion pa- 
per |28j an analysis of the quantum regime revealed, in 
the absence of dissipation, a C v ~ T 3 behaviour. 

An elastic system, such as the vortex lattice (VL) with 
external field aligned with z axis, is described by a N- 
component vector displacement field u a (Ri,z) (2V = 2 
for the VL in d = 3). The equilibrium positions Ri form 
a perfect TV-dimensional lattice of spacing a. Interactions 
result in an elastic energy H e i associated to the phonons 
of the vortex lattices H e i[u] = \ J q u a (q)^ a ^(q)up(—q). 

Here I g = f B z jkp f-t% W) denotes integration on the 
first Brillouin zone, q = (q±, q z ) a d dim- vector, and s is 
the distance between layers. For the triangular VL: 
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with P^Jk) = k a k p /k 2 and P T = I - P L . The disper- 
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sion of elastic moduli is implicit whenever needed. Impu- 
rity disorder is modeled by a short range gaussian random 
potential with in plane correlator A(r) = 8se^e~ r /( 4r /) 
interacting with the local vortex density. Here 2rf = £ 
the superconducting coherence length, e = (<I>o/47rA) 2 
is the vortex energy scale per unit length along z and 
8 a (small) dimensionless disorder parameter |35 j. The 
equilibrium Bragg glass phase (absence of dislocations, 
a/R a <C 1, R a being the translational correlation length) 
is described by the replicated partition function Z n = 
Tre~P HeS , (3 = l/T and — denotes disorder average. Af- 
ter standard manipulations 29] the replicated Hamilto- 
nian becomes H c g[u] = ^2 a H e i[u a ] + Hdi S [u] with: 



Ha 



I d 2 rdz R(u a (r, z) 



u b (r,z)) (2) 



Here R(u) = p\ J2k ^ k cos (K ■ u) in terms of po the av- 
erage vortex density and the disorder harmonics Ak = 



Id 



N ue iK-u 



A(u) at the reciprocal lattice vectors. More 



generally, an elastic manifold (such as a directed polymer 
d = 1) in a TV-dimensional embedding space in presence 
of a random potential W(u, x) is described by a simi- 
lar model with W(u, x)W(u', x) (x — x')R(u — it'), 
R(u) = —NV(u 2 /N). We compute the specific heat per 
unit volume C v (T) = — ^ where F is the free energy, 
£1 = L d the volume. We present the method using an 
isotropic disorder and elasticity tensor ^ a p(q) — cq 2 5 a p, 
and generalize to the vortex lattice later. 

To obtain the low T behavior, a first approach would 
be to assume a single minimum of the energy H [u] and 
expand around it. One then finds: 



C V (T) = C eq + AT + 0(T 2 ) 



(3) 



The exact expression for the linear term A = Ai mini 
given in |30| , involves cubic and quartic anharmonicity in 
H [u] in a given disorder realization. Disorder averaging is 
only easy to perform perturbatively in disorder, yielding: 



A = A 



lmin,pert 



jt(vtyR(u)\ u=0 



(4) 



with Ji = j q l/(cq 2 ) (for the manifold problem it gives 
A = AJl{N 2 + 6N + 8)V"'(0)/(3N)). Although such a 
single minimum low T expansion is useful for pure sys- 
tems, such as the Sine-Gordon model, for disordered sys- 
tems more than one minimum typically exists beyond the 
Larkin length R c . The resulting contribution to the spe- 
cific heat can then be estimated combining the droplet 
picture with the two levels argument 0, |^|. At each 
length scale / each of the subsystems i S [L/l) d may 
have a low lying secondary minimum (droplet) at excita- 
tion energy Ei independently distributed with probabil- 
ity P(E)dE = if (tO^IB!) where Ec is the typica 1 
pinning energy E c — cr^R^ 2 , 8 the free energy expo- 
nent. Approximating the contribution from scale / to the 



specific heat as: 
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treating as independent two level systems, C v = 
Ii>r T^i 1S dominated by the smallest scales, yielding: 
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While the two wells-droplet argument estimates as 
Adrop the contribution only of scales larger than R c , one 
can only hope to use A\ m i n to estimate the contributions 
of scales smaller than R c . It is thus instructive to com- 
pare them. The perturbative expression Q is infrared 
divergent for d < 2 as thermal fluctuations diverge. If 
one restricts by hand the integral in (@J to q > 1 /R c one 
finds a contribution of the same order than the droplet 
one Adrop- In d > 2 the integral is instead controlled 
by small scales and the droplet contribution is then sub- 
dominant. The two levels droplet model of Q can be 
improved by including anharmonicity in each well, for 
identical wells it simply adds A = A\ m i n + Adrop- 

The variational method [2^, extends these phe- 
nomenological considerations into a first principle quan- 
titative calculation in which the Larkin length naturally 
appears in a self consistent way. We introduce a Gaus- 
sian trial Hamiltonian H = \ f„G~p ab (q)u^(q)u b p(-q) 
which minimizes the variational free energy F val = Fq + 
(H c ff — H )h , where F denotes the free energy calcu- 
lated with Hq. The specific heat C' v = lim rw o 
can be reexpressed as: 

Cv = dT ^ (( E H <* W)He« + 2(H dis [u}) Hee ) (6) 

a 

Note the factor of 2 due to the (3 dependence of the dis- 
order term. Here we evaluate these averages using the 
variational hamiltonian Hq instead of the exact one H c s- 
Thanks to the variational equations this is equivalent to 
C v (T) = -(T/n)d 2 F var /dT 2 

We applied this variational approach to some pure 
models and checked that it is quite accurate [3(| . At 
low T it exactly matches the expansion around the mini- 
mum A var = Aimin and we checked, for the Sine-Gordon 
model in d = 2, that it is identical to the low T expansion 
of the exact result [33|. 

In the disordered case the solution of the variational 
equations requires Replica Symmetry Breaking (RSB). 
One denotes G(q) = G aa {q) and parametrizes G a ^b(ci) 
by G(q,u), where < u < 1 and similarly for B a b(x 
0) = B ab = (K(x) - 
For d > 2 one finds 
breakpoint u c and 



l2 )/iV with B = and B{u). 
a continuous RSB with a 



B(u >u c )=B = 2 7 TJ!(E) 



(7) 
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1 = -4 7 F"(S)J 2 (E) 
1 



Jn{z) 



(cq 2 + z) 1 



cr: 



(8) 
(9) 
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where 7 = 1, V(B) = A K exp(-6^ 2 /2) and 

R c is the Larkin length [36(. Eq. © is the so called 
marginality condition (MC) which also holds for the one 
step solution in d — 2. Starting from the expression: 



J q 



du[V(0)-V(B(u))] (10) 



which, as B(w) = B(u) (setting w = u/T) and w c = 
u c /T 37], turns out to depend implicitly on T only 
through £ and B, and using H7I9|) and ©, one obtains 
for the specific heat 



C v (T) = C eq +^F(B) 



(11) 



F{B) = V(B) - 7(0) - BV'(B) + -B 2 V"(B) 



, ,, - 2- ,v /" — J 9 is the equipartition value, 



where C eg = \M/Vt 

i.e. half the total number of modes per unit volume 
pjif . Eq. Hll|) is valid for any T, for periodic objects 
(Bragg glass) as well as manifolds, and independently 
of the replica structure of the solution (if broken, pro- 
vided MC holds). One thus finds that disorder increases 
the specific heat which has now a maximum and de- 
creases back to equipartition C eq at high T. Expand- 
ing (|11J) at low T one finds again © with an amplitude 
A var = § ! NV'"(0)J 1 (Y, T=0 ) 3 . For weak disorder R c > a, 
dispersionless elasticity and r/ < a, one finds from Ij7l9|l 
A var ps A/(cr 2 R c a 3 ), A — " ff ^ +4 ^ indeed larger by a 
factor (R c /a) 3 than 0. It also confirms the above dis- 
cussion: thanks to RSB (i.e. E / 0), the problems (e.g. 
in d = 2) of single minimum (i.e. replica symmetric) per- 
turbation theory are cured, the Larkin length being the 
natural scale. A plot of C V (T) is shown on Fig.^ 

The generalization to the vortex lattice using and 
N = 2 is straightforward. Eq. (f 1 1|) still holds, but in the 
formulaes (|7I8|I which determine B and £ = c§qR~ 2 one 
sets 7=5 and replaces J n (£) = J g (c669l + c 449z + £)~™ 
(we neglect compression modes, i.e. cee/cn 1). 

Taking dispersion in C44 and anisotropy into account 
yields several regimes analyzed in |30j , and for simplicity 
we present here only results for the weak disorder regime 
R c > a. The fluctuation of a vortex position is then mea- 
sured from B = TJi(S) w TJi(0) = c 2 L a 2 T/T m , where 
T m is the melting temperature and cl ~ 0.1 — 0.2 is the 
Lindemann number. Inserting in l|ll|l this yields C V (T) 
which has a maximum for T = T* i.e. B = B* determined 
by solving AF{B*) = {B*fV"'{B*). These fluctuations 
are estimated as Ji(S) ~ l/(4a v /c66C44) where here and 
below we denote C44 = 044(5^ = n /a) (below the dimen- 
sional crossover field, ^/{c^a 2 ) <C 1/s 2 , the integral 
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FIG. 1: Specific heat C V (T) — C eq in arbitrary units as a 
function of T/T* . 



over q z can be extended up to infinity and the integral is 
dominated by q± = it /a). 

We can now discuss in details the behaviour of C V {T) 
for rf < a. In this limit V(B) takes the dependence 
V(B) m -D/(2rj + B), with D = 5r 2 se 2 /a 2 , for B < a 2 



which holds up to melting. One finds: 

C V {T) - C eq + AT '/(l + T/(2T*)) 3 
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(12) 
(13) 



The maximum occurs for B* = r 2 and thus T — T* — 
Tmr 2 / (c 2 L a 2 ) , i.e. the so called depinning temperature 



, which can be below melting. The amplitude A (and 
large T tail) is independent of the detailed form of V(B) 
and is strongly enhanced by the anisotropy and disper- 
sion of C44. The value at the maximum is C* — C eq — 
8AT*/27. Using c 44 « e eVa 2 , c 66 w e /(4a 2 ) the 
calculation of J2(E) shows 30] that R c > a holds for 

8 < 6 C — 4-7T 3 / 2 — J. At this value of the disorder one 
obtains A » c 2 /V m r 2 ae) and C* - C e? » l/(a 3 e) ( 
e < 1 the anisotropy parameter [Jj). C* — C e? should 
be compared with the equipartition value, estimated as 
C eq = l/(sa 2 ). These become comparable around the 
dimensional crossover field B = B cr such that s ~ at. 

The above classical contribution l|12(l will hold only 
above the vorton Debye temperature , below which 
quantum effects become important. depends on the 
vortex mass, the Hall force and the friction force, which 
arise from the coupling of moving vortices with the nor- 
mal electrons bath in presence of scattering. Estimates 
for range j2^ from 10~ 3 T C in dirty (friction domi- 
nated) materials to lO" 1 !^ in supcrclcan limit (Hall dom- 
inated). There should thus exist a broad regime of tem- 
perature and field where the Bragg glass is stable and the 
result O holds. 

To assess whether this contribution l|12|) is observable, 
let us compare it with other terms linear in T in C v . First 
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the normal electrons in the vortex core lead to [l4j|: 



@v corc(-^) 



Tfk 



c2 



(14) 



where Tf is the Fermi temperature of the normal metal, 
a = 1 for s-wave superconductor |39j and a = 1/2 for 
lines of nodes in the gap |l5j. Given the large ratio 
Tf/T m = O(10 2 ) for e.g. YBaCuO, one finds, comparing 
A core and A<n s that the contribution from the cores can 
be comparable or smaller than the one from the disorder. 

There are other contributions from the vortex lattice 
as well. The mean field specific heat being C m f ~ 
e T/(T 2 a 2 ), the ratio A dls /A mf « ' {c 2 L aT c /{r f T m )) 2 
can be large. We have also computed the contribu- 
tion from non linear elasticity of the VL. Performing, 
as in |34l| . an expansion of the vortex interaction energy 
/ d> z YUj e oKo((Ri - Rj + Ui{z) ~ Uj(z))/X), we obtain: 



2 „2 



,CtQ 



A n i = Po^oa ' (-^— ) (74 + l3VcMC 66 ( 



T, 



T 

- 1 - rn. 



a 



(15) 



up to 0((a/X) 2 ) where the numerical prefactors 74 and 
73 are complicated lattice sums, given in ,30]. One can 
then estimate A n i/Adi S ~ r 2 /a 2 and this contribution 
is likely to remain small until melting. None of these 
contributions |4*oj is expected to exhibit a maximum at 
scale distinct from T m . 

To conclude, we found by explicit calculation that the 
vortex lattice classical contribution to C v due to dis- 
order can be important at the very least in the range 
10 _1 T C to T m , and possibly more. The coefficient Adi S 
of the low T linear behavior is magnetic field dependent, 
Adis oc l/(T m a). It is convenient to express results us- 
ing the meltin g, tem perature T m , which is experimentally 
measured |24j, |23, |2j| and can also be estimated from 
the standard elastic expression of T m = 4a 3 ^ 'cqqc^c 2 l 
allowing to extract directly the magnetic field depen- 
dence. T m ~ 1/VB leads to T* oc \f~B and Ad ls oc B. 
C V (T) exhibits, as compared to other contributions, a dis- 
tinct maximum around the dcpinning temperature scale, 
whenever smaller than melting. A crude estimate of T* 
is T* ~ T m B/(c 2 L H c2 ), which for a typical c L ~ 0.12 
gives a T* which is a fraction of T m for fields of about 
a Tesla (a fraction of order unity for YBaCuO around 
B ~ 10T the so-called tricritical point [13] )• It would be 
very interesting to perform precise measurements of C v 
to check the present proposal. 

We thank T. Klein, P. Monceau, A. Junot and CM. 
Varma for useful discussions. 
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